PASJ: Publ. Astron. Soc. Japan , 1-??, 
© 2013. Astronomical Society of Japan. 



Non-Thermal Radio and Gamma-Ray Emission from a Supernova 
Remnant by the Blast Wave Breaking Out of the Circumstellar Matter 

Takafumi Shimizu 1 , Kuniaki Masai 1 , and Katsuji Koyama 2,3 
1 Department of Physics, Tokyo Metropolitan University, 1-1 Mmami-Ohsawa, Hachioji, Tokyo 192-0397 

t-shimizu@phys. se. tmu. ac.jp 

2 Department of Physics, Graduate school of Science, Kyoto University, Oiwake-cho, Kitashirakawa, Kyoto 606-8502 
£C) ' 3 Department of Earth and Space Science, Graduate School of Science, Osaka University, 1-1 Machikaneyama, Toyonaka, 

Osaka 560-0043 

O 

' (Received ; accepted ) 

5—( ! 

, Abstract 

Wc calculate synchrotron radio emission and 7-ray emission due to bremsstrahlung, inverse-Compton 
scattering and 7r°-decay from the remnant of supernova which exploded in the circumstellar matter (CSM) 
formed by the progenitor's stellar wind. This sort of situation is a possible origin of mixed-morphology 
supernova remnants (SNRs) like W49B, which exhibit recombination-radiation spectra in X-ray emission. 
We assume that the CSM of 1.5 M Q exists at 0.07-3 pc away from the supernova in the interstellar medium 
(ISM) of density 0.016 cm" 3 . When the blast wave breaks out of the CSM into the ISM, its velocity rapidly 
increases and hence particle acceleration is enhanced. The maximum energy of protons reaches ~ 1300 TeV 
just after the break-out with ~ 0.5% of the explosion energy. Wc consider the non-thermal emission from 
the blast-shocked ISM shell after the break-out. Synchrotron radio flux at 1 GHz is tens Jy, comparable 
to the observed from mixed-morphology SNRs. Because of low density, the 7-ray luminosity is dominated 
by inverse-Compton scattering, which is higher than the 7r°-decay luminosity by an order of magnitude. 
» I The total 7-ray luminosity including bremsstrahlung is of the order of 10 33 ergs -1 lower than the typical 

value 10 35 -10 36 ergs -1 observed from mixed-morphology SNRs. However, if, e.g., ~ 10% of accelerated 
protons interact with some matter of density ~ 100 cm -3 , 7r°-decay 7-ray luminosity would be enhanced 
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to be comparable with the observed value. 
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^ ; et al. 2010d, Abdo et al. 2010c, Castro & Slane 2010). 

• . About 70% of Galactic supernova remnants (SNRs) ex- This may imply that such intense 7-rays come from the 

hibit shell-like morphology in radio wavelength (Green dense matter around MM SNRs. Actually, for these MM 

[ 2009). The radio shells are due to synchrotron emission SNRs, interaction with molecular clouds is suggested by 

■ by non-thermal electrons of order of GcV. Many of these 1720 MHz OH maser (Frail et al. 1994 , Yusef-Zadeh et al. 

• • I SNRs also exhibit shell-like morphology in X-rays, due to 1995, Frail et al. 1996, Green et al. 1997, Claussen et al. 

. ■ bremsstrahlung and line emission by thermal electrons of 1997, Hewitt & Yusef-Zadeh 2009) and/or near-infrared 

order of keV. The electrons are heated and/or accelerated observations (Keohane et al. 2007). The 7-ray spectra of 

$H . by interstellar shocks of SNRs. On the other hand, ~ 10% observed SNRs are not always single power-laws but cx- 

of the SNRs that have the radio shells exhibit center-filled hibit a break at an energy of ~l-5 GeV above which their 

thermal X-rays, and are called mixed-morphology (MM) slopes are steepen. 

SNRs (e.g., Rho & Petre 1998). Recently, Suzaku revealed X-ray characteristics different than shell-like SNRs and 

recombination radiation, which arc hardly expected for possible interaction with molecular clouds may suggest 

typical shell-like SNRs, in X-ray spectra of six MM SNRs that MM SNRs are remnants of core-collapse explosion 

so far observed: IC443 (Yamaguchi et al. 2009), W49B of massive stars surrounded by HII regions, stellar wind 

(Ozawaet al. 2009), G359.1-0.5 (Ohnishi et al. 2011), W28 matter, and molecular clouds, as in star-forming regions. 

(Sawada & Koyama 2012), W44 (Uchida et al. 2012), and Recombination radiation in X-rays are predicted by Itoh 

G346.6-0.2 (Yamauchi et al. 2012). & Masai (1989) for the remnant of a supernova sur- 

MM SNRs are characteristic also in 7-ray emission, rounded by its progenitor's stellar wind matter. They 

Fermi detected GeV 7-rays from SNRs: while the lumi- show the X-ray spectra due to rarefaction (adiabatic cool- 

nosities in the 1-100 GeV band of shell- like SNRs are ing) caused when the blast wave breaks out of the wind 

10 33 -10 35 erg s -1 (Abdo et al. 2010b, Castro & Slane matter to expand rapidly into the ambient interstellar 

2010, Katagiri et al. 2011, Abdo et al. 2011, Giordano medium. Recently, Shimizu, Masai & Koyama (2012), 

et al. 2011, Tanaka et al. 2011), the luminosities of MM hereafter Paper I, extend this work to non-spherically- 

SNRs arc distinctively higher, 10 35 -10 36 erg s -1 (Abdo et symmetric stellar wind matter, and find that rccombi- 
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nation X-rays exhibit center-filled morphology like MM 
SNRs with various shapes depending on the viewing di- 
rection. They also suggest that synchrotron radio shell is 
located outside, surrounding the X-ray emitting region. 

If rarefaction caused by the shock break-out is the origin 
of recombination X-rays found in MM SNRs, it is natu- 
rally of our interest whether the SNR model can explain 
observed radio and 7-ray emissions as well. Hence, in 
the present paper, we investigate non-thermal particles, 
which can be accelerated by the shock of the SNR model 
in Paper I, and emission thereby from the blast-shocked 
ISM shell. In the following section, we describe the SNR 
model in Paper I and the particle acceleration. In Section 
3 we show calculations of non-thermal radio flux and 7-ray 
luminosity, and discuss the results in Section 4. 

2. Model 

2.1. Supernova remnant 

We consider model B2 of Paper I for a model of MM 
SNR. The model describes evolution that the initially 
spherically-symmetric ejecta interact with anisotropic cir- 
cumstellar matter (CSM). Outside the CSM, a uniform in- 
terstellar medium (ISM) of density 0.016 cm -3 is assumed. 
Such low density can be possible for an H II region formed 
by the progenitor. The ejecta have an initial kinetic en- 
ergy of 2 x 10 51 erg and a mass of 10 Mq, and therefore 
the initial velocity of the blast wave is 8.5 x 10 3 km s" 1 . 
For comparison, we also calculate the evolution of SNR 
without CSM, which expands directly into the uniform 
ISM of the same density. 

The CSM is composed of the stellar wind matter blown 
by the progenitor in its pre-supernova phase. If the pro- 
genitor rotated rapidly, the stellar wind may have an 
anisotropic density distribution. We assume that the CSM 
is concentrated in the equatorial plane. The density on 
the equatorial plane is 4 times higher than that in the 
polar direction at the same radius. The mass-loss rate is 
5 x 10~ 5 M Q yr -1 at a wind velocity of 100 km s _1 . The 
inner and outer radii of the CSM are 0.07 pc and 3 pc in 
the equatorial direction. These radii imply that the wind 
activity lasts 3 x 10 4 years and then ceases 6 x 10 2 years 
before the explosion. The mass of the CSM is 1.5 Mq, 
which is obtained from the period of the wind activity 
and the mass-loss rate. 

The wind parameters can be possible for B[e] super- 
giants or luminous blue variables (LBVs). For exam- 
ple, radio observations suggest that W9, which is a B[e] 
supergiant in Westerlund 1, has the wind velocity of 
~ 10 2 km s _1 and the mass-loss rate of ~ 10 -4 Mq yr -1 
(Dougherty et al. 2010). LBVs have typically the wind 
velocity of ~ 10 2 km s _1 and the mass- loss rate of ~ 
10~ 5 - 10~ 4 M Q yr" 1 (Humphreys & Davidson 1994). 
Although LBVs were considered not to explode as su- 
pernova, recent observation shows that a progenitor of 
SN 2005gl is a LBV (Gal- Yam et al. 2007). Also, ob- 
servations of type Iln supernovae suggest the interaction 
between ejecta and dense CSM, which have the wind ve- 
locity of ~ 10 2 — 10 3 km s _1 and the mass-loss rate of 



- 10" 4 - 1 Mq yr" 1 (e.g., Kiewe et al. 2012). 

The blast wave breaks out of the CSM at ~ 450 yr in 
the equatorial direction, and then is rapidly accelerated. 
Rarefaction wave propagates inward from CSM-ISM con- 
tact interface. This causes rapid adiabatic expansion and 
thus cooling of the once-shocked CSM and ejecta, and re- 
sults in recombination-radiation X-rays. Since the ISM 
is rarefied enough to make a density difference of factor 
~ 10 at the CSM-ISM interface, rarefaction occurs for 
the mass loss rate 5 x 10 -5 M© yr -1 assumed here (cf. 
Moriya 2012). After that, blast wave propagates to form 
a shocked shell into the ISM, while the second reverse- 
shock propagates inward. Figure 1 shows a snap shot, 
density map of the shocked matter averaged over the line 
of sight, of the model SNR. The inner black lines represent 
the reverse-shocked ejecta which is bright in thermal X- 
rays, and the outer gray lines represent the blast-shocked 
ISM which is faint in X-rays but can be bright in radio. 
The late time evolution of the blast wave approaches that 
of the model without CSM, as shown by Itoh & Masai 
(1989). 

Figure 2 shows the blast- wave radius i?t>, velocity V s , 
mean temperatures (T p ) of protons and (T c ) of elec- 
trons, and mean number density (n) of the blast-shocked 
ISM. The proton and electron temperatures are calcu- 
lated from the temperatures kT Pt0 cx m Pi eK, 2 at the shock 
front, assuming the energy transport from protons to elec- 
trons through Coulomb collisions in the post-shock region 
(Masai 1994; see also Paper I), where m p and m e are the 
mass of proton and electron, respectively. 




Fig. 1. Contours of the number density of shocked matter 
averaged over line of sight in the equatorial plane of model 
B2 at 700 yr. Black lines represent those for reverse-shocked 
ejecta. Gray lines represent those for the shocked ISM. "RR 
X-rays" means recombination-radiation X-rays. 

2.2. Particle acceleration 

We consider diffusive shock acceleration by the blast 
wave: part of thermal particles of the shocked ISM are 
scattered across the shock by magnetic inhomogeneities 
and gain momentum. The energy spectrum of accelerated 
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particles is expressed in a form (Bell 1978) 

N(E) = K(E + mc 2 )(E 2 + 2mc 2 E)-^ +1)/2 
with 



(1) 



K 

for E\ 



e(n)( M -l)(i? i 2 nj + 2m C 2 i? inj )^- 1 )/ 2 



<E<E ma *. Here, £=( 7 -l)? 



(2) 

2 is the kinetic 
energy of the particle of mass m with the Lorentz factor 
7, i?i n j is the injection kinetic energy, -E max is the maxi- 
mum kinetic energy of accelerated particles, and £ is the 
injection efficiency The injection efficiency is defined as 
the ration of the number density of accelerated to thermal 
particles. 

The diffusive shock acceleration results in a single 
power-law energy spectrum, as described above. On the 
other hand, 7-ray observations of SNRs suggest that the 
energy spectrum of particles is not simply a single power- 
law, but, for instance, Fermi-observed SNRs show a break 
at an energy of ~l-5 GeV, as mentioned in Section 1. 
Since cooling time at the break energy is much longer 
than the age, this break may reflect acceleration processes. 
Therefore, for the energy spectrum of accelerated particles 
we assume a broken power-law 



N(E) = 



KE~», 
KE~^ 



l E~»\ 



for E < E h , 
for E > E h , 



(3) 



where Eh is the break energy of 10 GeV, taken so as to 
make a GeV break in the 7-ray spectrum. The spectral 
indexes are assumed to be \i = 2 and fi2 = 2.3, which is 
a medium value of spectral index of cosmic-ray sources 
(e.g., Putze et al. 2011). 

For the injection energy £"i n j and efficiency £, we con- 
sider that particles in the high energy tail of the thermal 
distribution are injected into acceleration process, as 

-^inj — p,c Ap C /cXp_ G . (4) 

with a constant A, where the characters "p" and "e" mean 
proton and electron, respectively. Then a relation between 
£ and A is given by 



JZ°'f p , (E)dE 



l-erf(A p , c 1/2 ) + 



1-1/2 



p.c 



l/2 -A p 



(5) 



where f(E) is the Maxwellian distribution function, and 
erf is the error function. 

We determine £ p for the pressure of accelerated protons 
to be equal to 10% of the ram pressure of ISM enters in 
the blast wave. The pressure of accelerated particles is 
given by 



Pi 



CR 



Pinj 



:£(n)cpinj 



In 



N'{p)pvdp 
2pb 



rnc 



PbC 



-M2+2' 



(6) 



where we use /i = 2 in the last expression. Here, v is 
the particle velocity, p is the momentum of a particle, 
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Fig. 2. From top to bottom, evolutions of the blast-wave ra- 
dius, the blast-wave velocity, the proton and electron tem- 
peratures averaged over the blast-shocked ISM in the equa- 
torial direction, and the number density averaged over the 
blast-shocked ISM in the equatorial direction, as functions of 
the elapsed time after explosion. 

Pinj = (2m£ in j) 1/2 is the injection momentum, 
E max /c is the maximum momentum, p^, is the break mo- 
mentum, and N'(p)dp = N(E)dE . In the last expres- 
sion in equation (6), pi n j <C mc and p max ^ Pb > rnc 
arc considered. The injection efficiency of protons is 
roughly proportional to the blast- wave velocity V s because 
£ p oc V s 2 /p in j„p oc V s 2 /Tp /2 oc V s . The injection efficiency of 
protons reaches the maximum ~2x 10~ 4 at~530 yr, and 
then decreases to 5 x 10~ 5 at ~ 10000 yr. We determine 
£ e for the pressure -PcR-e of accelerated electrons not to 
exceed the pressure -Pcr-p of accelerated protons. The 
ratio of the pressure of accelerated electrons to protons is 

1/2 

(7) 



-PcR-p 



0.05^ 



E- ■ 

Einj-p 



If the injection energy of electrons is the same as protons, 
£e ^5 20 £ p follows. In the following, wc express £ c in unit 

The maximum energy -E max is determined by the time- 
scales of energy gain and loss. Adiabatic loss due to SNR 
expansion is negligible through the age concerned here. 
The dominant loss process is synchrotron radiation and 
inverse- Compton scattering for electrons. Assuming that 
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1) mean free path of a particle is its gyration radius (Bohm 
limit), 2) shock compression ratio is 4, and 3) accelerated 
particles are relativistic (7 3> 1), we estimate the time- 
scales of acceleration and radiation loss as 

t a .^^£ (8) 



SeBV 2 



and 



Qnm c 



^loss(electron) — / D p rr \j (9) 

respectively, where <tt is the Thomson scattering cross 
section, e is the elementary electric charge, B is the 
strength of the magnetic field in the shock downstream, 
assumed to be 4 times the strength in the upstream, and 
C^cmb is the energy density of cosmic microwave back- 
ground. In SNRs, the magnetic field strength can be 
stronger than the average interstellar value by magnetic 
amplification mechanisms, as suggested by X-ray variabil- 
ity of RX J1713. 7-3946 (Uchiyama et al. 2007). Using the 
equation (14) of Bell & Lucek (2001), we calculate the field 
strength in the SNR evolution. The magnetic amplifica- 
tion may cause the non-linear feedback from accelerated 
particles to the shock structure. However, such feedback 
is small when the injection efficiency is lower than ~ 10~ 4 
(e.g., Fcrrand et al. 2010), which is marginally attained 
after the blast-wave break-out concerned here. 
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Fig. 3. Maximum energy that the accelerated protons and 
electrons can reaches at a given time (upper) and the magnetic 
fields strength in the shock downstream (lower), as functions 
of the elapsed time after explosion. Gray line represents the 
maximum energy of the protons in the SNR evolution without 
CSM. 

In Figure 3 we show the time evolution of E max and B. 
For protons E max is determined by t acc ~ t age and reaches 



- 1300 TeV at - 700 yr, while - 800 TeV at - 430 yr in 
the case without CSM. For electrons E max is determined 
by t acc ~ iioss(eiectron) j and its maximum is about 10 TeV 
at the moment of the break-out. At ~ 700 yr just af- 
ter the break-out, E max takes its maximum/minimum for 
protons/electrons because of rapid increase of the shock 
velocity and magnetic field. For the explosion energy 
2 x 10 51 erg assumed, the total energy of accelerated pro- 



tons is 1 x 10 49 erg at 
- 10000 yr. 



700 yr and - 2 x 10 5U erg at 



3. Non-thermal radiation 

We calculate the non-thermal emission from the blast- 
shocked ISM shell, in which accelerated particles are con- 
fined for the either shorter time, t agc or R 2 /D. Here D 
is the diffusion coefficient and taken to be in the form 
D W (E/10 GeV)(B/10 (J,G)- X cm 2 s" 1 with a numerical 
coefficient Dio. Observations of cosmic-rays suggest that 
D~ 10 28 cm 2 s" 1 at 10 GcV (Berezinskii et al. 1990). On 
the other hand, near SNRs, GcV and TeV observations 
suggest D ~ 10 26 cm 2 s" 1 at 10 GeV (e.g., Li & Chen 
2012). We adopt O 10 = 3x 10 27 so that D ~ 10 28 cm 2 s" 1 
for B = 3 /iG, typical field in the interstellar space, and 
D ~ 10 26 cm 2 s _1 for B ~ 100 /iG, which could be attained 
for SNRs. 

3.1. Synchrotron radio 

Synchrotron radiation at the frequency v = 1 GHz 
is mainly emitted by electrons with energies ~ 
2 (v/l GHz) 1 / 2 (5/100 A*G)- 1/2 GeV. According to 
Ginzburg & Syrovatskii (1965), a flux F syn (v) of syn- 
chrotron radiation emitted from electrons with the broken 
power-law spectrum is 



F f 



syn(^) 

1 



2m r 



3e 



-(M-l)/2 



F{x)x^ )/2 dx + E-^ 



3eB 



x / F(x)x^ 2 - 3)/2 dx 



Airr 2 dr 1 



(10) 



where 



16m^c 5 ^ 
X= 3eBE 2 ' 

K e is K of electrons (see Eq. (3)), d is the distance to the 
SNR, r is the radius from the center of the SNR, and F{x) 
is the synchrotron function. The characters "inj" , "b" , 
and "max" are correspond to the injection, break, and 
maximum energy, respectively. Integration interval of r is 
given by the shocked ISM shell defined in the beginning 
of Section 3. Figure 4 shows the time evolution of the 
1 GHz flux for the magnetic field in Figure 3. Since the 
blast wave is little decelerated, the radio flux continues to 
increase through ~ 10000 yr (see Section 4). 



Non-thermal emission from SNR with CSM 
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Fig. 4. Left: Radio flux at 1 GHz of synchrotron radiation from the blast-shocked ISM shell as functions of the elapsed time after 
explosion, calculated for £ c = 10 £ p and d = 8 kpc. Right: Luminosities of brcmsstrahlung (broken), invcrsc-Compton scattering 
(dotted) and 7r°-decay (solid) 7-rays from the blast-shocked ISM shell in the 1-100 GeV band, as functions of the elapsed time after 
the explosion, calculated for £ = 10 £ p . 



3.2. Bremsstrahlung j-ray 



3.3. Inverse-Compton ^-ray 



Relativistic electrons emit bremsstrahlung 7-rays by 
interacting with target protons. Number of the 
brcmsstrahlung photons emitted from electrons with the 
broken power-law spectrum per unit time per unit energy 
per unit volume is (Blumenthal & Gould 1970) 



dN~ 



dtdhvdV 
Aar 2 cn^ 



hv 



In 



Ahv 



dEN c E- 



-E 2 
3 



4 

- -Ehv- 
3 



4ar cK c tit 





1" 






\ m e c z J 


~ 2 



(3^ 2 +n+i){hvy 
3m(m-1)(//+1) 

i^-^E^ _ ( M2 



(hvf 



+ 



CC! 1 

(3M2+M2+4)(fei/)-^ 



1)0*3 + 1) 

hv > ±E h , 



-i) 
hv < 



;E h 



(11) 



3m 2 (M2-1)(M2 + 1) ' " 2 J 

where a is the fine structure constant, r is the classical 
electron radius, h is the Planck constant, and nx is the 
number density of target protons. The target protons are 
assumed to be thermal protons in the shocked ISM shell. 
We consider E max 3> hv in the last expression in equation 
(11). Using equation (11), we get a bremsstrahlung 7-ray 
luminosity by 



Lh\-c 



hv- 



dN~, 



dtdhvdV 



-47rr drdhv. 



(12) 



Integration interval of r is given by the shocked ISM shell 
defined in the beginning of Section 3. We calculate the 
luminosity using (n 2 ) obtained from hydrodynamical cal- 
culation of Paper I. Figure 4 shows the time evolution of 
7-ray luminosity due to bremsstrahlung in the 1-100 GeV 
band. 



When a photon of energy eo is scattered by a relativistic 
electron of Lorentz factor 7, scattered photon has energy 
of ~ 7 2 eo- Because the maximum energy of accelerated 
electrons is ~ 10 TcV, the electrons can emit photons of 
energies up to ^100 GeV by scattering cosmic microwave 
background (CMB) photons. Number of photons emitted 
from electrons with the broken power-law spectrum by 
scattering CMB photons of the temperature Tcmb per 
unit time per unit energy per unit volume is (Blumenthal 
& Gould 1970) 



dN~, 



dtdhvdV 

_ 87r 2 r 2 K c E-^ 



h 3 c 2 (m e c 2 )^ 2 
2^+ 3 ( M 2 + 4/i 2 



— (fcXcMB) 
ID 



(M2+5)/2 



(hv 



-( M 2 + l)/2, 



(/i2 + 3) 2 (/i 2 + 5)(/i 2 + l) 



fJ-2 



M2 + 5 



-hv- 



2^+ 2 7T 2 



3(^2 + 1) \AE^ kTcun 



(M2+l)/2 



(13) 



where £ is the zeta function and E maXie is the maximum 
energy of accelerated electrons. Using equation (13), we 
get a inverse-Compton 7-ray luminosity by 



Lie — 



hv- 



dN~, 



dtdhvdV 



-Airr drdhv. 



(14) 



Integration interval of r is given by the shocked ISM shell 
defined in the beginning of Section 3. Figure 4 shows 
the time evolution of the 7-ray luminosity due to inverse- 
Compton scattering in the 1-100 GeV band. 

S.^. tt° -decay j-ray 

Relativistic protons emit neutral 7r°s through inelastic 
collisions with protons, and then the 7r°s decays into two 
7-ray photons. Wc calculate the 7r°-decay 7-ray lumi- 
nosity, using the parameterized cross section of inelastic 
proton-proton collision 
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[Vol. 



ci(-Ep) — 3 



0.95 + 0.06 In 



xl0~ 26 cm 2 , (15) 



i Gcvy_ 

and (5-function approximation of number of 7r°s emitted 
per unit time per unit energy per unit volume 



dtdE^dV 

CUT 



■CTinol m„C 




(16) 



which are used in Aharonian & Atoyan (2000). Here 
f n ~ 0.17 is mean fraction of the kinetic energy of pro- 
ton transferred to ir° per collision, E p = r ym p c 2 is total 
energy of proton and E„ = jm^c 2 is total energy of 7r°. 
Number of 7r°-decay photons emitted from protons with 
the broken power-law spectrum per unit time per unit 
energy per unit volume is 



dN~ 



dtdhvdV 

= 2 
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{El - m 2 .c 4 ) 1 /2 dtdE^dV 

1 *"« -*L 



-dE, 
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(17) 



for hv < fnEh, 



(£)' 



M2 

for hv > f„Eb, 



0.95 + 0.06 In 



1 GoV 



where A m i n = hv + (to 2 c 4 /4/i^) is the minimum pion en- 
ergy to produce photon of energy hv and K p is A' of pro- 
tons (see Eq. (3)). In the second expression in equa- 
tion (17), hv > m n c 2 and E^ > m^c 2 are considered be- 
cause we calculate photons above 1 GeV. In the last ex- 
pression of equation (17), we approximate the proton en- 
ergy spectrum as the relativistic form and the variable 
of the spectrum as (m p c 2 +E w /f w ) ~ E^/f^, because 
En/fir > (1 GcV/0.17) - 6 GeV > m p c 2 . Using equation 
(17), we get a 7r°-decay 7-ray luminosity by 



hv- 



dN. 
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dtdhvdV 



-Airr drdhv. 



(18) 



Integration interval of r is given by the shocked ISM shell 
defined in the beginning of Section 3. As in the calcula- 
tion of the brcmsstrahlung luminosity, we use (n 2 ) instead 
of (n) 2 . Figure 4 shows the time evolution of 7-ray lumi- 
nosity due to 7r°-decay in the 1-100 GeV band. 

4. Discussion 

For the low density ISM of density 0.016 cm -3 in the 
present model, supposed for an HII region (e.g., formed 



by the progenitor and extended to a few tens pc), the blast 
wave is little decelerated through ~ 10000 yr. As a result, 
in the context of diffusive shock acceleration described in 
Section 2, the radio flux continues to increase, because the 
increase of the emission measure overcomes the decrease 
of the magnetic field strength. Consequently, for about 
ten thousand years, recombination-radiation X-rays are 
observed from the irregular-shape inner part of SNR (see 
Figure 1 and Paper I) , while the radio emission of tens Jy 
is observed from the blast-shocked ISM shell. 

In the beginning of the Sedov/Taylor phase where the 
blast wave is being decelerated significantly as V s oc i -3 /' 5 , 
the radio flux turns to decrease slowly as i~ 3 / 10 , and 
then approaches nearly constant as the magnetic field ap- 
proaches its interstellar value (~ 3 /xG) and T e approaches 
T p . Also the inverse- Compton 7-rays turns to decrease as 
cx t -1 / 5 in the Sedov/Taylor phase, while 7r°-decay 7-rays 
are nearly constant. This sort of analysis is done also for 
the phase < 10000 yr with the relation V s oc t~ s where s 
is given by the hydrodynamical calculation, and gives a 
good agreement with the computed time evolution of the 
radio and 7-ray emission in Figure 4. It should be noted 
that s ~ 0.4 at 10000 yr, yet smaller than the Sedov value 
s = 3/5, and the SNR is in the transient phase to the 
Sedov/Taylor regime. 

Again because of low density, the 7-ray luminos- 
ity of the shocked ISM shell is dominated by inverse- 
Compton scattering through the SNR evolution con- 
cerned. However, 7r°-dccay 7-rays could be enhanced by 
interactions with dense external matter, e.g., dense HI 
gas, molecular clouds or a cavity wall formed by the stel- 
lar wind of the progenitor. If 10% of accelerated protons 
interact with such matter of density n ~ 100 cm -3 , the 
luminosity L T would exceed 10 35 erg s _1 at a few thou- 
sands year, comparable to the typical 7-ray luminosity of 
MM SNRs. The interactions with molecular clouds are 
suggested in many MM SNRs by OH maser and/or near- 
infrared observations. The interaction with HI gas is sug- 
gested in RX J1713. 7-3946 by observations (Fukui et al. 
2012), and may be expected also in MM SNRs. 

Interaction with some dense external matter may be 
realized on the 7-ray to radio flux ratio. We show the ratio 
of 1-100 GeV to 1 GHz flux in Figure 5. One can see that 
the ratio is systematically high for MM SNRs compared to 
shell-like SNRs except for RX J1713. 7-3946 and Vela Jr. 
For RX J1713. 7-3946 interaction with molecular clouds 
(e.g., Dame et al. 2001) and/or HI gas (Fukui et al. 
2012) is suggested by observations. As for Vela Jr., an 
X-ray source CXOU J085201.4-461753, possible neutron 
star, is located near the center of the SNR (Pavlov et al. 
2001). That might be related to the hard radio spectrum 
and high 7-ray/radio ratio, though the pulsar activity is 
not observed. The low ratio of Cas A is due likely to a 
strong field ~ 1 mG (e.g., Arbutina et al. 2012). 

When electrons responsible for radio and elec- 
trons/protons for 7-rays are produced in the same volume, 
the 7-ray/radio ratio is reduced to ut/B 3 / 2 . Therefore, 
if the particle spectrum and B do not vary so much from 
SNR to SNR, the ratio can be a measure of the density 
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Fig. 5. Flux ratio of 1-100 GcV to 1 GHz of shell-like (open 
circles) and MM (filled circles) SNRs observed so far by 
Fermi. Attached numbers indicate the power indexes of the 
radio spectra (Green 2009) and "RR" means that recombi- 
nation-radiation X-rays are observed. 

References— Radio flux of SNRs, except for RX J1713.7-3946 
(Acero et al. 2009); Green 2009, 7-rays and distance: Tycho's 
SNR; Hayato et al. 2010, Giordano et al. 2011, G349. 7+0.2; 
Frail et al. 1996, Castro & Slane 2010, CTB 37A; Reynoso & 
Mangum 2000, Castro & Slane 2010, Cygnus Loop; Blair et 
al. 2005, Katagiri et al. 2011, Cassiopeia A; Reed et al. 1995, 
Abdo et al. 2010b, Vela Jr.; Katsuda et al. 2008, Tanaka et 
al. 2011, RX J1713.7-3946; Fukui et al. 2003, Abdo et al. 
2011, IC443; Welsh & Sallmen 2003, Abdo et al. 2010c, 
W49B; Moffett & Reynolds 1994, Abdo et al. 2010e, W28; 
Velazquez et al. 2002, Abdo et al. 2010d, W51C; Koo et al. 
1995, Abdo et al. 2009, W44; Wolszczan et al. 1991, Abdo 
et al. 2010a, W30; Fich et al. 1989, Castro & Slane 2010, 
3C 391; Frail et al. 1996, Castro & Slane 2010 

n-T of the matter with which the particles interact. For 
B ~ 100 nG the high ratios observed from MM SNRs may 
be explained by 7r°-decay if the density of the target mat- 
ter tit > 10 cm~ 3 , higher than the typical ISM density 
1 cm -3 ). The high ratios could be explained also by 
inverse-Compton if B < 10 fiG. Such a field may be possi- 
ble for shell-like SNRs, but unlikely for MM SNRs which 
exhibit rather high radio flux. 

Finally, we mention the effect of the CSM, stellar wind 
matter here. An important effect of the CSM is that 
the shock break-out raises the maximum energy E^x to 
~ 1300 TeV for protons (see Figure 3). Since E niax oc 
BV s 2 t cx V s 3 oc (75 oj /i\/ cj ) 3 / 2 , where 75 oj and M oj are the 
initial kinetic energy of ejecta and the ejecta mass, re- 
spectively, E max would reach ~ 3000 TeV, the cosmic-ray 
knee energy, for 2 times larger value of E e j/M e j than that 
in the present model. 

We are grateful to Yutaka Ohira, Inouc Tsuyoshi and 
Ryo Yamazaki for meaningful discussion about particle 
acceleration. Also to the anonymous referee for his/her 
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Research 22540253 and 24540229, from Japan Society for 
the Promotion of Science (JSPS). 
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